N92-2193 8 

COMPTEL Solar Flare Observations 

J.M. Ryan 3 , H. Aarts 2 , K. Bennett 4 , H. Debrunner 5 , C. de Vries 2 , J.W. den Herder 2 , G. 
Eymann 4 , D.J. Forrest 3 , R. Diehl 1 , W. Hermsen 2 , R. Kippen 3 , L. Kuiper 2 , 

J. Lockwood 3 , M. Loomis 3 , G. Lichti 1 , J. Macri 3 , M. McConnell 3 , D. Morris 3 , V. 
Schonfelder 1 , G. Simpson 3 , M. Snelling 4 , H. Steinle 1 , A. Strong 1 , B.N. 
Swanenburg 2 , W.R. Webber 3 , and C. Winkler 4 

^ax-Planck Institut fur Extraterrestrische Physik, Garching, Germany 

2 Laboratory for Space Research Leiden, PB 9504, 2300 RA Leiden, The Netherlands 

3 Space Science Center, Institute for the Study of Earth, Oceans and Space University of New 
Hampshire Durham, NH 

4 Space Science Division of ESA, ESTEC, Noordwijk, The Netherlands 
^University of Bern, Bern, Switzerland 

COMPTEL as part of a solar Target of Opportunity campaign observed the sun during the 
period of high solar activity from 7 June to 15 June 1991. Major flares were observed on June 
9 and 1 1. Although, both flares were large GOES events (> X10), they were not extraordinary 
in terms of y-ray emission. Only the decay phase of the June 15 flare was observed by 
COMPTEL. We report the preliminary analysis of data from these flares, including the First 
spectroscopic measurement of solar flare neutrons. The deuterium formation line at 2.223 MeV 
was present in both events, and for at least the June 9 event, was comparable to the flux in the 
nuclear line region of 4 - 8 MeV, consistent with SMM observations. A clear neutron signal 
was present in the flare of 9 June with the spectrum extending up to 80 MeV and consistent in 
time with the emission of y-rays, confirming the utility of COMPTEL in measuring the solar 
neutron flux at low energies. The neutron flux below 100 MeV appears to be lower than that of 
the 1982 June 3 flare by more than an order of magnitude. The neutron signal of the 1 1 June 
event is under study. Severe dead time effects resulting from the intense thermal X-rays require 
significant corrections to the measured flux which increase the magnitude of the associated 
systematic uncertainties. 

1 . Introduction 

As a secondary scientific objective of the COMPTEL y-ray telescope on the Arthur 
Holly Compton Gamma Ray Observatory, observing and measuring solar flares assumed an 
unexpected importance due to the unusually high level of solar activity early in the mission. 
During the two week period of 1991 June 1 to June 15 region 6659 (Boulder reference number) 
transited the disk producing an unprecedented number of large flares. After the first week of 
this activity, the prospect of more flares was such that the sun was declared a Target of 
Opportunity and on June 7 the observation of the Cygnus region was suspended and the 
Compton Observatory was reoriented so as to place the sun near the spacecraft z axis. Since the 
spacecraft z axis coincides with the axis of the COMPTEL field-of-view, this enabled 
COMPTEL to utilize its full capabilities in measuring the gamma ray and neutron fluxes from 
the flares which occurred the following week. Major flares were recorded on June 9, 1 1 and 15 
after which the Compton Observatory re-assumed its original objective of the all-sky survey. 
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2 COMPTEL v-Ray and Neutron Response 

COMPTEL is a two element Y -ray <^rrw^ 
used to measure the photon s en ^g y andinc Qn scattering? neutron detection relies 

measure neutrons. However, ^ther thanusi g ; C Pj radiation are important in studying high 

SrXXc p“ I' descrip, ion of ,he insmrmem ard its in-flight performance ,s 
given by den Herder et al. (1991). 



COMPTEL 

IMAGING COMPTON TELESCOPE 

Fifrnre 1 Schematic of COMPTEL showing the interactions of 
Y-rayT and neutrons. The scat.er angle ♦ for Y-rays ,s grven by* 

= cos-Vl - e/£2 +e/(E\ + E 2 )', where e = 511 keV. The non 
relativistic hard sphere, neutron-proton scattering formula 

tan^ <p = E\IE S . 

The basic scattering process for photons (and ^even element forward 

quantities that COMPTEL the locatio n of the scatter in the 

SIS n ttS&S tte energy deposit E 2 in that detector. One computes 
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from these interaction locations the direction of the scattered photon's (neutron's) velocity 
vector. For y-rays, one computes the scatter angle 0and the total y-ray energy (E\ + E 2 ), 
where 


0 = cos - ] (1 - eJE 2 + {/(El + E 2 )). 

Here, eis the electron rest mass energy and <p is the Compton scatter angle provided E\ + E 2 is 
the full incident y-ray energy. Without measuring the direction of the scattered electron in Dl, 
only the polar angle and not the azimuth angle of the scatter is known. The time-of- flight (TOF) 
is also measured and used to help identify y-ray scatters; whereas for neutron detection, it is 
used to measure the energy of the scattered neutron. 

Much like ay-ray event, the ideal type of neutron interaction in COMPTEL occurs when 
the incoming neutron elastically scatters off a hydrogen nucleus in the Dl detector. The 
scattered neutron then proceeds to the D2 detector where it interacts, depositing some of its 
energy to produce a trigger signal as indicated in Figure 1. The energy of the incident neutron is 
computed by summing the proton recoil energy Ei in the Dl detector with the energy of the 
scattered neutron E s deduced from the TOF from the Dl to the D2 detector. The scatter angle 
for non-relativistic neutrons (< 100 MeV) can be computed by the formula: 

tan 2 <p = E\!E S . 

Relativistic corrections become important above 100 MeV, in which case the expression 
becomes 


tan 2 (l> = 


EllEs 

1 + En/2m n ' 


where E n = E\ + E s is the incident neutron kinetic energy. 

The material in D 1 is a liquid organic scintillator, NE213A, with the properties of low 
density and low Z (H/C ratio = 1.286). The material elastically scatters both y-rays and 
neutrons off atomic electrons and hydrogen nuclei, respectively. The detector thickness has 
been chosen such that it is a fraction of a mean free path thick, meaning that the incident y-ray 
(or neutron) can scatter in Dl and usually leave Dl without scattering again. 

The liquid organic scintillator in Dl (NE213A) possesses pulse shape discrimination 
properties, in that energetic protons produce light pulses with longer rise times than those of 
electrons (and other minimum ionizing particles). This capability allows for efficient 
identification of recoil proton signals produced by fast neutrons elastically scattering off 
hydrogen in Dl. Any reaction producing a recoil proton can be identified by these means, such 
as inelastic scattering of fast neutrons off carbon producing a y-ray and either a knock-on proton 
or neutron (which can then elastically scatter off hydrogen). Pure y-ray producing neutron- 
carbon reactions in Dl also occur and represent an intrinsic background in both y-ray and 
neutron measurements. 

Since only the scatter angle 0 is measured, the incident particle direction is constrained to 
a cone mantle of half angle 0. For y-rays, errors in the measured energy and placement of the 
cone mantle occur via uncertainties in the measured energies in Dl and D2 and uncertainties in 
the measured interaction positions. Partial energy absorption in D2 (from an escaping y-ray) 
yields a low value for the total y-ray energy and a large value for the scatter angle 0. Similar 
errors are inherent in the neutron measurement as well, except that the TOF resolution is the 
primary component in the scattered neutron energy uncertainty. The neutron energy resolution 
(a) of about 6% results in an uncertainty in the production time of the neutron at the sun of - 1 
minute for a 40 MeV neutron. 
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Figure 2. 

The energy spectrum from a monoenergetic y-ray source 
produced under the constraint that the y-ray scatter angle * be 
consistent with the source direction. 

Similarly in detecting and measuring neutrons from the sun the neutrons can be traced 
, f tr, ni thmiiuh the ancle <t> onto a cone mantle restricting the incident 

Sec" to include the sum This 

requiring a scatter angle <p which is consistent with the solar direction. 

4. y-ray Flare Overview 

y-ray emission was detected from at least three Hares observed by ( te^ntenshy 
June 199 1 solar Target of Opportunity obtKwatton. ^own m F.gure 3 bJ» we j, £ lhe 

resolution to show the 
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apparent motion of the sun over the course of two days. Figure 4 is the COMPTEL image of 
the sun on June 9 during the impulsive phase of the flare. The angular resolution of COMPTEL 
is on the order of one degree while the diameter of the sun is one half that, so the images do not 
reveal any information on the location of the y-ray emission on the solar disk. The flare on June 
15 was not observable until about 0900 UT, approximately 40 minutes after the X-ray 
maximum. Before this the spacecraft was in orbit night or in the South Atlantic Anomaly. 
However, an image of the sun could still be constructed, supporting the conclusion of Akimov 
et al. (1991) of prolonged (> 1 hour) y-ray emission from this flare. 




Figure 3. The intensity-time profile > 1 MeV and the time- 
integrated count spectrum from the solar flare of 1991 June 9. 
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The intense solar flare on June 4 was also measured by COMPTEL, although the sun 
^ r- * ac' TUt* rnMPTFT hurst detector system measured the y ray 

^r^ 

ofjreutrons which was detectable a. ground .eve, a, 
Mt. Norikura in Japan (Muraki et cil. 1991). 

5 . Analysis of the 1991 June 9 Solar Flare 

Only the June 9 flare which occurred during the Target of Opportunity period has been 
. 1 1 Ttif* cmoivck includes both Y-ravs and neutrons. 
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Figure 4. The COMPTEL image of the sun taken during the 
impulsive phase of the flare. The coordinates are declination as 
Right Ascension. The image is consistent wjthapointsource 
response. The correct solar position is a - 5h 07m and d 
22.9°. 

5a. y-Ray Analysis 

TU T Qfl^nn Yinoi occurred in active region 6659 at N34E04 with the X-ray 
The June 9 flare (2 , * j ^ nT p t u e COMPTEL data we define the impulsive 

star, a, 0134 UT and ^maximum a, 0143 UT. From he ^COMEIb show „ in Figure 3 

phase to be from 0136 to 0146 UT Y y & V d £ t which has a threshold 
These data were obtained from the COMF1 bLjngn range ou d 5ackground COU nt 

of ~ 1 MeV The broad maximum near 0158 UT is due only to tn e increaseu uac s 

sills 
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the impulsive phase the y-ray count rate was high, filling up the yl telemetry channel and 
precluding the transmission of any y2 events. During this same time interval the activity in the 
forward veto detector VI was high as well. This was due to the intense thermal X-ray flux at a 
few keV piling up in VI, continuously triggering that detector and producing a large and 
variable dead time effect in COMPTEL. The thermal X-ray flux remained high, well beyond the 
impulsive phase. This required additional analysis to estimate the dead time for both the 
impulsive phase and the period of time afterward during which we expect the arrival of solar 
neutrons. The dead time correction was estimated by the throughput of LED calibration events 
in D1 during the impulsive phase and the period afterward. During the impulsive phase the 
average live time fraction was (28 ± 1)%, while from 0146 to 0205 UT the effect was (78 ± 
4)%. 

In the analysis of the y-ray data it is important to estimate both the integrated flux of the 
2.223 MeV emission from deuterium formation and the integrated flux in the nuclear line region 
of 4 - 8 MeV. The emission of the 2.2 MeV line depends upon the number of free, low energy 
neutrons within the photosphere, w'hile emission in the broad nuclear line region depends upon 
the energy spectrum of protons above 30 MeV as they interact with the solar atmosphere. Since 
the production cross sections of neutrons and nuclear lines have very different energy 
dependences, the comparison of the y-ray emission at 2.2 MeV and in the range of 4 - 8 MeV is 
a useful indicator of the parent proton spectral shape (Murphy and Ramaty 1984). 

The broad band efficiency from 4 - 8 MeV was estimated by a Monte Carlo simulation 
calculation under the same restrictions as imposed on the real data. The instrument efficiency in 
the line at 2.2 MeV was obtained through interpolation of calibration simulations, again with 
restrictions similar to those applied to the real data. Using both the live time estimates and the 
estimates of the telescope efficiencies, we calculate from the time-integrated count rate spectrum 
(Fig. 3) the time-integrated fluxes to be 32 and 35 y-cnr^ for the the y-ray emission in the range 
of 4 - 8 MeV and in the 2.223 MeV line, respectively. The uncertainties for these quantities are 
estimated to be 50 and 30%, respectively. The uncertainties are largely determined by the 
systematic errors in the estimated efficiencies rather than by dead time corrections or statistics. 
No other spectral analysis has been performed, although there is evidence for other identifiable 
lines in the spectrum. 

5b. Neutron Analysis 

Additional processing must take place in the neutron data analysis. All the quantities 
normally measured and telemetered are in terms of y-ray interactions. Quantities appropriate to 
neutrons must be computed. These include (1) converting the electron-equivalent energy 
measured in Dl, the upper liquid organic scintillator, into the proton-equivalent energy, (2) 
calculating the energy of the scattered neutron (from Dl to D2) from the time-of-flight, (3) 
computing the scatter angle </>from the n-p scattering kinematics, (4) discarding events which 
have electron-recoil pulse shape signatures in Dl and (5) discarding events which fall within the 
y-ray time-of-flight window. The neutron differential flux can then be calculated based upon 
the data satisfying these criteria along with the live time and the neutron detection efficiency. 

Presently, the best figures for the COMPTEL neutron detection efficiency are the 
extrapolations of the COMPTEL prototype neutron measurements to the full COMPTEL 
configuration. The calibrations were performed at only a few energies and incident angles 
(Ryan eta/. 1987). The efficiencies needed for the flux calculations were interpolated between 
calibration points. The differential flux was calculated taking into account the limited number of 
operating detectors and the live time during the neutron measurements. 

In principle the data restrictions listed above coupled with data restrictions that are 
usually used for y-ray spectral studies should be all that is required to derive a solar neutron 
spectrum and intensity-time profile. The main selection criterion needed for spectral studies is 
that the scatter angle 0 of an event be consistent with the position of the sun. This means that 
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the event circle whose center is defined by the velocity vector of 'the .scatter, fri neutron and ^ 
whose diameter is defined by the scatter a "J ’ = c'orrect'si nee in Compton telescopes the imaging 
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consfderable and could 11 Wgh as a factor o/two. Better figures for these correcttons, m 

for the loss of neutrons through f decay over the 1 AU flight path. The resu hug dtfeen tal 
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from the measurement of neutron-decay protons in interplanetary space. Based upon the 

difference between our value and the spectrum of Evenson et al., we can conclude that there 

should be no measurable signal in any ground level neutron monitor as there was for the 1982 

June 3 flare, which was measurable at the few a level. 

6. Conclusions 

(1) The sun was imaged in y-rays for the first time. 

(2) The first direct energy measurement of solar neutrons was performed. 

(3) The y-ray fluence in the nuclear line range of 4-8 MeV was approximately 30 y-cnr 2 . 
Within errors this equals the fluence in the 2.2 MeV line, making it consistent with many 
flares measured by SMM (Forrest, priv. comm.). 

(4) The integrated neutron spectrum from 15 to 80 MeV is on the order of 2 x 10 28 n-sr 1 
at the sun. 

(5) The neutron emission was consistent in time with the y-ray emission, i.e. from 0136 - 
0150 UT. There is strong evidence against precursor neutron emission. Prolonged 
emission after the impulsive phase appears to be consistent with background. 

(6) The neutron background measured by COMPTEL was surprising large and required 
severe data cuts to eliminate, reducing the effective area accordingly. Consequently, we 
have revised our estimate of the neutron sensitivity to ~ 6 n-cnr 2 for a 20 minute 
observation, or equivalently, 5 x 10' 3 n-cm^-s' 1 in the range of 15 to 80 MeV. 
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